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Abstract 



We study the production of isolated photons in e + e~ annihilation and give 
the proof of the all-order factorization of the collinear singularities. These 
singularities are absorbed in the standard fragmentation functions of partons 
into a photon, while the effects of the isolation are consistently included in the 
short- distance cross section. We compute this cross section at order as and 
show that it contains large double logarithms of the isolation parameters. We 
explain the physical origin of these logarithms and discuss the possibility to 
resum them to all orders in as- 
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1 Introduction 



At LEP, the Tevatron and the future LHC, the detection of high-energy photons produced 
by short-distance interactions goes through the definition of an isolation criterion that aims 
at decreasing fake photon signals coming from ir° decays. The high-energy tt°, which is 
seen as a single cluster in a calorimeter, belongs to a jet and is accompanied by hadronic 
energy, whereas a bremsstrahlung photon emitted by a quark is unaccompanied as soon as 
the photon-quark angle is large enough. Therefore any criterion to select 'unaccompanied 
photons' enhances the signal-background ratio. 

The criterion that is mostly used in current practice is the following. In e + e~ collisions 
at LEP, a photon of energy E 1 is said to be isolated if it is accompanied by less than a 
specific amount Ef l E 1 of hadronic energy E^ one in a cone of half-angle 5 around the photon 
momentum. This isolation criterion can be written as follows: 

E™* = y £E i Q(5-e* Y )<e h K r , e h >0, (1.1) 

where the energies and the relative angles 9^ are defined in the centre-of-mass frame. 
In hadron collisions at the Tevatron, the isolation criterion is similar, but the energies 
are replaced by the transverse energies while the angles and the cone are defined in the 
azimuth-pseudorapidity plane, so that 6 iry is replaced by = ^f]^ + 0? r Typical values 
of the isolation parameters and 5 are = (9(10 _1 ) and 5^0.7. 

The criterion in Eq. ( |1 . 1|) can be stated in an equivalent way by saying that the fraction 
of electromagnetic energy inside the isolation cone has to be larger than a fixed value x c : 



E. 



7 



>x c , (1.2) 



X^ 7 Ei Q(5 — 6> i7 ) + E.y 
where 

x c ee — L- < 1 . (1.3) 
1 + e h 

In the case when no isolation criterion is applied, the inclusive photon cross section is 
computable by using the QCD factorization formula: 



a 
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J2 C A (^ 2 ) ® D J/A (^) , (1.4) 



where (g> denotes the convolution over the energy fraction z. The hard-subprocess cross sec- 
tion Ca(^ 2 ', z), which describes the production of a high-energy parton A, can be calculated 
order by order in QCD perturbation theory. The non-perturbative phenomena are included 
in D 7 / a (z,/i 2 ), the inclusive fragmentation function of the QCD parton a (a = q,q,g) into 
a photon Jl|, , and the direct term (A = 7) not proportional to the fragmentation function 
is taken into account through the convention D 1 / 1 {z, fi 2 ) = 5(1 



z 



In the case of an isolated photon, the isolation criterion (11. If) enforces additional phase- 
space restrictions. This implies that the cross section is no longer fully inclusive and, hence, 
that the factorized expression (|1.4j) is not necessarily valid. 
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The applicability of the factorization theorem to isolated photons is thus a basic issue 
that has to be dealt with before the corresponding cross sections can be studied within the 
conventional QCD framework. The perturbative calculations that have been performed so 
far do not help to this purpose. In the case of hadron collisions, all the available calcula- 
tions || beyond the leading order (LO) in the strong coupling as are based on approximate 
methods that, in particular, lead to an incomplete treatment of the fragmentation contribu- 
tions. The next-to-leading order (NLO) calculation of Ref. H for e + e~ collisions introduces 
by definition a fragmentation component that explicitly depends on the isolation param- 
eters and thus differs from the process-independent fragmentation function in Eq. (|1.4j) . 
Actually, in the case of e + e~ annihilation, factorization has recently been questioned by 
Berger, Guo and Qiu ||. The results of Ref. || have been criticized by some of us in 
Ref. §. 

In the present paper we confirm that factorization is fulfilled at NLO and extend this 
argument by proving the validity of factorization to all orders in perturbation theory. 

Although we show that isolation does not spoil factorization, yet the phase-space re- 
strictions due to the isolation criterion ( [[ . 1\ ) are not harmless, as discussed in detail in 
the second part of the paper. The factorization theorem deals with collinear singularities 
that occur in the calculation at the parton level. Once these singularities, whose origin is 
non-perturbative, have been absorbed in the fragmentation functions, the short-distance 
cross section can still have a divergent behaviour at some points of the phase space when 
computed order by order in perturbation theory 0. These divergences are due to certain 
kinematical constraints that, limiting the fully- inclusive character of the cross section, pro- 
duce an imperfect compensation between real and virtual emission of soft (and collinear) 
partons. In isolated-photon cross sections, the soft-gluon divergences are double logarith- 
mic and appear at a specific point inside the phase space || || . Owing to its perturbative 
origin, this disease can be cured by summing the logarithmic divergences to all orders in 
perturbation theory. 

The outline of the paper is as follows. In Sect. ^ we give our proof of the validity 
of the factorization theorem for isolated-photon cross sections defined by the criterion^ in 



Eq. We show how the effects of the isolation are consistently included in the short- 

distance subprocess and, in particular, we discuss the functional dependence of the short- 
distance cross section on the isolation parameters. In Sect. |3] we consider isolated photons 
produced in e + e~ annihilation and we compute the NLO contribution to the fragmentation 
component of the short-distance cross section. We present results in analytic form for any 
value of the isolation parameters Eh, 8. Using these explicit expressions, in Sect. § we discuss 
in detail the physical origin of the soft-gluon divergent behaviour of the NLO cross section 
in the vicinity of point x 1 = x c . Finally, in Sect. [| we summarize our results and outline 
how an all-order resummation can eventually lead to well-behaved theoretical predictions 
for the cross section. 



§The case of jets containing isolated photons has been considered in Refs. (J, |[ ||. An alternative 
definition of the isolated photon has recently been suggested by Frixione JhJ] . 



2 



2 Factorization 



The factorization issue regards both hadron and photon distributions. Thus, we make 
no distinction between these two cases, although isolated-hadron cross sections are of less 
experimental interest. 



2.1 e + e annihilation 

The customary factorization formula for the inclusive distribution of a single particle H 
with four-momentum p 1 produced in e + e~ annihilation is: 

1 da(Q 2 , x~) y-^ f 1 dx 



cr dx 
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W - D H/A (x 1 /x,fi 2 )Cf in) (a s ( f i 2 ),Q 2 /^x) + 0((l/Qr) , (2.1) 

A J X-y X 



where the sum extends over A = q/,q/,g when the observed particle if is a hadron and over 
A = qf, qf, g, 7 when H = 7 is a photon (in this case D 1 / 1 (z, jj 2 ) = 5(1 — z), by definition). 
The term 0({1/Q) P ) on the right-hand side denotes corrections that are suppressed by 
some power p > 1 of the centre-of-mass energy Q when Q 3> Aqcd- 



The formula (|2.1| ) states that at large values of Q 2 and for any fixed value of the energy 
fraction x 7 = 2p 7 • Q/Q 2 , all the long-distance physics phenomena can be absorbed in 
the non-perturbative fragmentation function D h /a(z, h 2 ) of the parton A into the particle 
H. The remaining coefficient function C^ ull \as(^ 2 ), Q 2 / fi 2 ; x) is short-distance-dominated 
and depends only on the partonic subprocess. 

The predictivity of Eq. ( |2.1[ ) within perturbation theory follows from the fact that not 
only C}± but also the (Revolution of the fragmentation functions are perturbatively 
computable as power-series expansions in a s . After having extracted D h /a(z, Ql) from 
experimental data at a certain value Q 2 = Qq, perturbative QCD predicts the cross section 
in Eq. ( |2.1| ) for any other value of Q 2 . 



The content of the factorization theorem in perturbation theory is nonetheless wider. 
It states that similar factorization formulae are valid for other less inclusive observables 
and that, in these formulae, the dominant non-perturbative contribution is universal and 
accounted for by the same fragmentation functions D h /a(z, /i 2 ) as in Eq. (|2~T|). Thus, the 
factorization formulae are obtained from ( |2.1|) by the replacement 

Cf ll \a s ^ 2 ),Q 2 /^x) - C A (a s (fi 2 ),Q 2 /^x,{J}) , (2.2) 

where {J} denotes the dependence on the particular observable. The differences among 
the various observables only regard the perturbatively computable coefficient functions Ca 
and the size and the power p of the power-suppressed corrections. 

A particular class of observables to which the factorization theorem applies is the class 
formed by what we call jet-type observables. These observables can be easily identified by 
examining how they are defined (and measured) in terms of the momenta of the final-state 
particles in the process. The definition has to fulfil the requirements of i) infrared safety, 
ii) collinear safety, and Hi) collinear factorizability. 
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The first two requirements regard the dependence on all particle momenta but the trig- 
gered momentum p 7 : infrared safety means that the value of the observable is independent 
of the momenta of arbitrarily soft particles, and collinear safety implies that, when some 
final-state particles are produced collinearly, the value of the observable depends on their 
total momentum rather than on the momentum of each of them. Collinear factorizability 
means that from the measurement of the observable one cannot distinguish whether the 
triggered momentum p 7 is carried by the particle H or by that particle accompanied by a 
bunch of particles parallel to it. 



These properties can be stated in a formal way as follows ITT], [12], [L3|]. Let us first 
introduce the exclusive cross section da!^ xcl \p^,pi, . . . ,p n ) to produce n + 1 particles with 
momenta p~,,pi, . . . ,p n , so that the single-particle inclusive distribution in Eq. (|2.1|) can be 
written in the following form 



da(Q 2 , 



= W da^ cl \p„ Pl , ...,p n ) 5{ Xl - 2 Pj ■ Q/Q 2 ) , (2.3) 

UXj n JQ(p~ f ,pi,...,p„) 

where the integration extends over the full (n + l)-particle phase space f2(p 7 , »i, . . . ,p n ). 
According to this notation, any less inclusive cross section daj is given by 



daj(Q 2 , x 



E / da^ cl -\p 7 , Pl ,...,p n ) 



dx~f n JQ(p 1 ,pi,...,p„) 

■ S(x, - 2 Pj ■ Q/Q 2 ) Fi n) (Q, Pl , {J};p a , ... Pn ), (2.4) 

where, for any exclusive final state with n + 1 particles of momenta Pj,pi, . . .p n , we have 
denoted by 

Fj n \Q,P~,,{J};Pu- ■ -Pn) (2.5) 
the measurement function that defines the actual observable. 

In terms of this function, the properties of a jet-type observable are£] 

i) infrared safety: 

Fj n+1) (Q,p 7 , {J}; p i ,...p n+1 ) -> Fj n) (Q,p 7 ,{J};pi, ...p n+x ) , (2.6) 



Pi -» o 



ii) collinear safety: 

Fj n+1) (Q,Pj, {J}; Pi, ■ -PuPj, ■ ■ -Pn+i) -* Fj n) (Q,p 7 , { J};pi, ...,pi + pj,.. .pn+i) 

Pi II Pj 



(2.7) 



in) collinear factorizability: 
Fj n+ ] {Q,p ,{J};p 1 ,...p i ,...p n+1 ) Fj n) (Q,p 7 +Pi,{J};Pi, Pn+i) ■ (2.1 



Pi ii p-, 



^In Eqs. ( |2.6| )-( ^8[ ) bold-face characters are used just to emphasize the differences between left-hand 
and right-hand sides. 
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In the case of the isolated-particle distribution, the measurement function is specified by 
the isolation criterion in Eq. ( [Lip . One sees that the expression ( |1 . 1| ) fulfils the properties 
( p.6|) and (|2.7| ), provided the isolation parameter is kept finiteQ. 



The fulfilment of the collinear factorizability in Eq. ( |2.8| ) may appear more problematic 
by naive inspection of Eq. ( |1 . 1[) . If the isolation parameters J = {eh, 5} were assumed to 
be the relevant variables for the short- distance subprocess, the measurement function could 
be defined as 



F( n+1 \Q, Pl , J = {e h , 5}; Pl , . . . Pn+1 ) = 6 



n+l 



1=1 



(2.9) 



Then, considering the collinear limit where, for instance, p n +i is parallel to p 7 : 



F^ +1 \Q, P7 ,J = {e h ,5}; Pl ,... Pn+1 ) 







thE-f — E n+ i —^Ei 0(5 — 6. 



\ 



F (n) {Q, Pl - Pn +x/e h , J = {e h , 5}; Pl , ... Pn ). (2.10) 



Since the function F^> on the right-hand side of Eq. ( |2.10|) depends on the momentum 
P-y—Pn+i/th rather than onp 7 +p n+ i, Eq. ( |2.8p is not fulfilled by the definition in Eq. (|2.9|). 



This effect can easily be understood. In its collinear-fragmentation process the triggered 
parton (particle) turns out to be necessarily less isolated (its energy decreases while the 
amount of accompanying hadronic energy increases) and, eventually, the isolation criterion 
can be violated. Thus one cannot insist on factorizing a short-distance subprocess that 
depends on the fixed isolation parameter e^. Hard partons produced at short distances have 
to be more isolated than the triggered particle H and their isolation has to be increased as 
the energy of H decreases. 

This is the key point to show that the isolated cross section is a jet-type observable. The 



isolation criterion in Eq. ( p. . 1|) can indeed be recast in a form that explicitly fulfils collinear 
factorizability by considering J = {r 7 , 5} as the relevant parameters of the hard-scattering 
subprocess. Here r 7 is defined by 



r 7 = — > 



2[j2E i e(5-9 i7 )+E 7 



Q 



> X. 



7 ' 



(2.11) 



and represents the upper limit on the total-energy fraction inside the isolation cone. Using 
these isolation parameters, the measurement function that corresponds to the criterion in 



Eq. (1J-) can be written as 



F^(Q,p„J = {r 1 ,5}; Pl ,... Pn ) = e 



Q 



^2 " 



V 



J2E i Q{5-e ij )+E 1 



1=1 



(2.12) 



"The violation of infrared safety in the case of perfect isolation, i.e. when E 1 ™ 110 = 0, was pointed out 
inRefs. 
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Equation Q2.8Q , as well as Eqs. ( |2.6|) and (|2.7|) , straightforwardly apply to the function in 



Eq. (|2.12|) as long as r 7 > 2Ey/Q = x T Note that, performing the limits of Eqs. (|2.6|) (|2T 



the parameters J = {r 7 , 5} in Eq. (|2.12|) have to be kept fixed and regarded as variables that 



are independent of the momenta Py,pi, . . -p n - In particular, this implies that the factorized 
coefficient functions of Eq. (|2.2|) explicitly depend on the parameters J = {r 7 , 5}. 

We thus conclude that, for the production cross section of isolated photons, factorization 
is valid and the all-order factorization formula is: 

1 da(Q 2 , x y , x c , 5) t-^ f 1 dx 



J2 I — D j/a (xy/x, ii 2 ) C a (a s (ti 2 ),Q 2 /iJ?;x,ry, 5) 

n - n r r, *> X>v X 



°"0 dXy a=q f ,qf,g "^7 

+ Cy(a s (n 2 ),Q 2 /n 2 ;Xy,ry,5) , (2.13) 

where C a and C 7 are the coefficient functions for the fragmentation and direct components, 
respectively. Note that the fragmentation function D 7 / a is independent of the isolation 
parameters 5 and Eh- The effect of the isolation is entirely included in C a and C 7 through 
their dependence on 5 and r 7 = Xyj (1 + Eh). 

A final remark is in order. At present, a field-theory proof of the factorization theorem 
is known only for fully inclusive deep inelastic scattering (DIS). In the case of less inclusive 
observables, factorization can be justified order by order in perturbation theory on the 
basis of general power-counting arguments ||15j| . The requirements of infrared and collinear 



safety and collinear factorizability are sufficient to guarantee the validity of factorization 
through power counting. In this respect, our proof of factorization for isolated-photon cross 
sections cannot be considered rigorous from a field-theory viewpoint, but it is certainly at 
the same level of rigour as for any other semi-inclusive cross sections [I5|. 



2.2 Hadroproduction and photoproduction collisions 

The discussion of the previous subsection can straightforwardly be extended to the produc- 
tion of isolated photons in collisions of hadrons and/or real photons. In these cases, the 
calculation of the cross section at the parton level contains additional singularities that are 
produced by initial-state collinear radiation. These singularities have to be absorbed in the 
non-perturbative parton distributions of the colliding particles. Thus, the essential differ- 
ence with respect to e + e~ annihilation is that, in hadronic collisions, the photon-isolation 
criterion must not spoil the factorization of the initial-state collinear singularities. 

As mentioned in Sect. [I], in the current experimental practice [[RJ the hadronic-collision 
versions of the isolation criterion (|1 . 1|) involve transverse energies and angular distances 
evaluated in the azimuth-pseudorapidity plane. These variables are invariant under lon- 
gitudinal boosts along the beam direction and, hence, they are insensitive to initial-state 
collinear radiation. It follows that these types of isolation criteria fulfil the QCD factoriza- 
tion theorem. 

Note, however, that the factorization of initial-state collinear singularities can be easily 
violated when boost-non-invariant variables are used. For instance, this is the case if one 
considers the angular distances _Rj 7 = yj r/f^ + 4>fy and insist in using energies rather than 
transverse energies. 
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3 Fixed-order calculation in e + e annihilation 



To explicitly check the factorization formula fl2.13|) for e + e~ annihilation, one needs a 
perturbative calculation to (at least) the first non-trivial order in as- We have carried out 
this calculation in analytic form and the results are presented in this section. We limit 
ourselves to the case in which the isolation parameters Eh, 5 vary in the range of practical 
interest, namely: 

< e h < 1 , < 5 < tt/2 . (3.1) 



The coefficient functions C a for the fragmentation component of Eq. ( |2.13| ) have the 
following perturbative expansions: 



+ 0(a 2 s ) , 
(3.2) 

where \ qf = \q f = e 2 f , \ g = Y,f e q f > C Qf = C qf = C q and e qf is the electric charge of the 
quark of flavour /. 

At the LO only the quark coefficient function contributes: 

Cf 0) (x,r 7 ) =6{l-x) 6(r 7 -l) , Cf 0) (x, r 7 ) = . (3.3) 

Note also that C^ LO \x,r^) is non-vanishing only for r 7 > 1, i.e. for a; 7 > x c . 
The results of our calculation for the NLO expressions can be written as 



Ci NLO \Q 2 /^x,r„8) = 9(r 7 



c mo,w ){QV ^, x) 



+ Ci NL0 ^\x,r,,6)] + Ci NL0 ^\x,r,,6) 



(3.4) 



where C^ ArLO,full - ) ((5 2 / /i 2 ; x) denote the customary NLO coefficient functions [|T7| , [T~8|] for the 
single-particle distribution, i.e. for the case in which no isolation is applied. Their explicit 
expressions in the MS factorization scheme are: 



Cf LO < fuI1 ) (Q 2 //i 2 ; 



x) 



1 +x 2 , (1 -x)Q 2 3 1 
In 



1 — x 
2 
3' 



\--^-^-)5{l-x) 



4 J 



2 1 -x 
5 3 
V 2 ~ 2 



+ 



1 + x 2 
1 — x 



In x 



x 



Cl NLO > Ml \Q 2 / l x 2 ;x)=C F 2 



l + (l-x) 2 (l-x)x 2 g 2 



In 



x 



(3.5) 



(3.6) 



The other two terms on the right-hand side of Eq. ( |3.4| ) derive from the isolation criterion 
and do not depend on the factorization scheme. We find that: 



The contributions C^ NLO ' m ^ are non- vanishing only in the region 

2 8 

r-y < 1 + tan - , 



(3.7) 
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and are given by 



c (NLO,in)fa rrr: S )=- Q(x+ - x)Q(x - X~) Ci 



1 + x 2 



X 



111 



(1 - r 7 + x)(l -x)tan 2 (5/2) 



- (4 - x) 



(3.1 



x — r 7 
1 — x 



+ 



x sin 2 (5/2) _ 



Cf LO ' in )(x,r 7 ,5) = -e(x+-x)e(x - x~) C F 
( 1 + (1 - xf (1 - r 7 + x)(l - x) tan 2 (5/2) 



(3.9) 



In 



x 



where 



r 7 — 1 



x — r~ + 



(1-x) 



1 -xsin 2 (5/2) 



4(r n 



^ r2 sin 2 (5/2) 
The contributions C^ NLO ' out ) are non- vanishing only in the region 



r 7 < 1 



and have the following explicit expressions 



C { q NLO ' out) (x,r^5) = 9(r 7 - x) 



1 + x 2 
1 — x 



In 



4 - x 



1 - x) tan 2 (5/2) 



1 l + xsin 2 (5/2) 
1 - x ~ 1 -x sin 2 (5/2) 



(3.10) 



(3.11) 



(3.12) 



Cf LO ' ont \x, r 7 , 5) = 9(r 7 - x) C t 
f 1 + (1 -x) 2 , 1 



In 



x 



x) tan 2 (5/2) 



- 2 



1 - 



-x)(l + xsin 2 (5/2)) 
1 - x sin 2 (5/2) 



(3.13) 



The origin of the various NLO contributions in Eq. (|3.4j ) is easily understood. We 
are interested in the process 7* — > q + q + g + 7 when one of the QCD partons is 
collinear to the photon. Owing to this collinear decay, to compute the coefficient function 
C( NLO \Q 2 1 /1 2 ; x, r 7 , 5), we simply have to evaluate the cross section^ for the three-parton 
subprocess 7* — > q + q + g when the triggered parton a (a = q, q or g) carries momentum 
p, parallel to p 7 , and energy fraction x = 2E/Q. We denote by Pi,P2 the momenta of the 
other two partons and by xi,X2 their energy fractions. According to the isolation criterion 
specified by Eq. (|2.12| ), the corresponding measurement functions is: 



F^(Q,p,{r J ,5}; Pl , P2 ) = e r 7 - £ x { 9(5 - 6 i7 ) 



x 



(3.14) 



i=l,2 



To make explicit the factorization of collinear singularities, we rewrite Eq. ( |3.14|) by adding 
and subtracting a contribution that is independent of the momenta, as follows: 



6(r 7 - 1) + [F (2 \Q,p, {r 7 ,S};pi,p 2 ) - 6(r 7 - 1) 



(3.15) 



'More technical details can be found in Ref. 
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When inserted into Eq. Q2.4|) and combined with the virtual correction, the first term in 
Eq. |U5|) gives exactly (cf. Eq. Q) the fully inclusive contribution Cf LO > full )(Q7/i 2 ; x) 
to Eq. 



Then, we have to consider the term in the square bracket in Eq. ( 3.15|) , which, on the 
basis of the general factorization argument of Sect. [2|, is expected to give a non-singular 
contribution. To show that, we decompose this term in two parts that correspond to the 
cases in which one additional parton, either p\ or p 2 , is inside the isolation cone (F^ 2,m )) 
and both partons are outside it (F( 2 > out )). Using Eq. ( gJ§ , we obtain: 



F®(Q,p, {r 7 ,6}- Pl ,p 2 )-Q(r 



= F^ in \Q,p,{r 1 ,6}- Pl ,p 2 ) 

+ F^ t \Q, P} {r 1 ,6};p l ,p 2 ) , (3.16) 



where 



F^(Q lP , {r 7 , 6};p u p 2 ) = - {&(6 - 0i 7 )[e(r 7 - 1) - 0(r 7 - x x - x)\ + (1 <-> 2)} 



F^(Q,p, {r 7 , 6}; Pl ,p 2 ) = 0(0 l7 - S) Q(6 2j - 5) [e(r 7 - x) - 6(r 7 - 1) 



(3.17) 
(3.18) 



Let us first consider the emission inside the cone. Since Xi+x > 1 because of kinematics, 
the term in the square bracket on the right-hand side of Eq. ( |3.17| ) corresponds to the phase- 
space region 

xi+x>r 7 >l. (3.19) 

This forbids the parton pi to become either soft {x\ = 0) or collinear to the photon (6*i 7 = 0) 
because in both cases the three-parton kinematics implies x\ + x = 1, thus violating the 
constraint ( p,19| ). Therefore, we can safely perform the integration over pi and obtain the 
finite contribution C[ NLO >' in \x, r 7 , 5) in Eqs. Q, ([D|), (|]|). Note that parton radiation 
inside the isolation cone is included in both terms C^ NLO ' ful1 ) and C^ NLO ' m ^ of Eq. ( |3.4j ). 
The contribution of C^ NLO ' m ^ is negative (cf. Eq. ( |3.17| )) and represents the suppression 
effect of the non-isolated distribution produced by the isolation criterion. 

A similar discussion applies to the emission outside the cone. Since x < 1 because of 
kinematics, the term in the square bracket on the right-hand side of Eq. Q3.18p vanishes 
when r 7 > 1. Thus the term in the square bracket can be replaced by 



1 > r 7 > x 



(3.20) 



In this region, p\ and p 2 cannot become either soft or collinear because in both cases 
one has x = 1. Therefore the integration over p\ and p 2 is safe and we obtain the finite 
contribution C^ LO ' out )(x, r 7 , 5) in Eqs. Q, flST^), (P7TB| ). Note that the C^ L °^ is 
positive, although smaller (cf. the subtraction in Eq. (|3.18|) ) than the full non-isolated 
contribution C^ NLO ' in the region r 7 < 1. 

We remind the reader that the coefficient function for the direct component of the 
factorization formula ( |2.13| ) was analytically computed to the lowest order by Kunszt and 
Trocsanyi M. Using our notation, their result can be written as follows: 



C 7 (as(/i ),Q In ;x 7 ,r 7 ,5) = —A 

Z7T 



— Cf LO \Q 2 /fi 2 ; x 7 , r 7 , 6) + 0(a s ^ 2 )) 



(3.21) 
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where a is the fine structure constant and C^ NLO \Q 2 / fi 2 ; x 7 , r 7 , 5) is given in Eq. 
The relation ( p.21|) between C 7 and our result for the NLO gluon coefficient function C^ NL °^ 



can be regarded as a partial check of the calculation described in this section. 



Our NLO results in Eqs. ( |3.4|) - (|3.13|) do not fully confirm those in Ref. ||. There, the 



coefficient functions C^ NLO > were found to be affected by singularities (single and double 
poles in 1/e, where e = 4 — d parametrizes the number d of space-time dimensions in di- 
mensional regularization) that would spoil conventional factorization at the specific point 
x 7 = x c . The method of calculation used in Ref. || and the interpretation of these singu- 
larites have been criticized in Ref. ||. The method described in this section clearly exhibits 
the factorization of the collinear singularities. Using the decomposition in Eq. ( |3.15|) , we 
separate a term (that in the square bracket), which is manifestly free from singularities^] 
and computable in 4 space-time dimensions, from a remaining contribution whose depen- 
dence on the isolation parameters is only due to the constraint r 7 > 1. This is an overall 
constraint (i.e. it does not act on the partonic variables) and its effect is thus harmless: the 
ensuing (dimensionally regularized) singularities are the collinear singularities that univer- 
sally enter the fully-inclusive cross section. 



4 Divergent behaviour for 



The NLO coefficient functions C { a NLO) (Q 2 / fi 2 ; x,r^, S) in Eq. (£§ are well-behaved for any 



value r 7 7^ 1, 



i.e. both for x 7 < x c and for x 7 > x c . However, when r 7 — > 1 they become 



divergent. The divergent behaviour (Fig. |T]) at this point x 7 = x c , which we shall call the 
critical point, can be easily derived from the explicit espressions in Eqs. ( |3.<j| ), Q3.9p and 



When the critical point is approached from below: 



(4.1) 



we obtain 



LO) (Q 2 //i 2 ; x, r 7 , 6) = 5(1 - x)C F [in 2 ((1 - r 7 ) tan 2 (5/2)) 



+ - ln((l-r 7 )tan 2 (<5/2) 



C^ LO \Q 2 /^x,r„5)=0(l). 



0(1) 



(4.2) 
(4.3) 



When the critical point is approached from above: 



x 



c — ) 



(4.4) 



^By 'manifestly free from singularities' we mean the following. If this term is evaluated in d = 4 — e 
space-time dimensions, then the limit e — > can safely be performed because it leads to the contributions 
Ci NLO ' ln ' (x, r 7 , 5), C'a NL °' ° ut ' (x, r 7 , 5) that are integrable in any interval of the kinematical variables 
x,r 7 ,<5. This integrability guarantees that the perturbative coefficient functions are well-defined distribu- 
tions. However, it does not imply that they are non-singular functions (cf. the discussion in Sect. |^). 
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X 



X. 



Figure 1: The divergent behaviour of the isolated-photon cross section in the vicinity of 
the critical point x 7 = x c at LO (solid), NLO (dotted) and as expected after resummation 
(dot-dashed). 




3 In' T 



2 Vtan 2 (5/2) / 

(4.5) 

l J \ mil yu i 6) j j 

C^ \Q^^,x,r 7 ,S) = C F 2 1 + {1 ~ X)2 In f -^-1 + 0(1) , (4.6) 



where 0(1) stands for any term that is finite for r 7 = 1. Owing to relation Q3.21 ), the 
coefficient function C 7 (as(/^ 2 ), Q 2 /fJ, 2 ; x 7 , r 7 , 5) for the direct component diverges when x c 
is approached from above. 

The presence of these logarithmic divergences was first pointed out in Ref. |J. It is not 
in contradiction with the factorization of the collinear singularities ||. Indeed, these two 
types of singularities have different physical origins, as discussed in rather general terms in 
Ref. J7J and recalled below. 

The factorization theorem deals with infrared singularities that affect the parton-level 
calculation for the whole range (or for finite intervals) of the relevant kinematical variables, 
e.g. x 7 , 5. These singularities are due to the long-distance component of the scattering 
process that is intrinsically non-perturbative, in the sense that it is not suppressed by some 
inverse power of Q as the hard-scattering scale Q increases. The validity of the factoriza- 
tion theorem for the isolated-photon cross section guarantees that this non-perturbative 
component can be described by the universal fragmentation function D 1 / a (x 1 Q 2 ). The 
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remaining contributions to the cross section are short-distance-dominated. They consist 
of i) power-suppressed terms (cf. the term 0((1/Q) P ) on the right-hand side of Eq. ( |2.1| )) 
that are controlled by non-perturbative phenomena and of it) a short- distance component, 
the coefficient functions Ca, which depends only logarithmically on Q and is computable 
as a power series expansion in as (Q 2 )- 

However, the fact that the coefficient functions are perturbatively computable does not 
imply that the coefficients of their perturbative expansion are non-singular functions. In 
fact, in general these coefficients are singular generalized functions or distributions that lead 
to finite quantities only when they are integrated with sufficiently smooth test functions. 
The 'plus'-distribution 

1 , 1 



w(x) = 2C7 



In 



1 — x 1 — x 



(4.7) 



that enters in Eq. ( |3.5| ) is a well-known example of this type of singular generalized func- 
tions. An analogous divergent behaviour, namely 



Cf LO \Q 2 /^, x, r 7 , 5) ~ e{l-rjeh-x)C F In T ^j T y j £ = ^' ( U 



is observed in the NLO coefficient functions of Eq. (|3.4j ) when the cone size S — > at fixed 
7 . The double- and single-logarithmic divergences in Eqs. (|4.2[ ), fl4.6p, ( |4.5|) are integrable 



r 



in any neighbourhood of the point x 7 = x c and essentially belong to the same type of 
singularities |TS[. These singularities are known as divergences of the Sudakov type and, in 



spite of the validity of the factorization theorem, they are still produced by the radiation 
of soft and/or collinear partons. 



The main difference between the divergent behaviour in Eqs. (|3.5|) and (|4.8| ) and that in 
Eqs. ( |4.2j ), ( |4.5| ), ( |4.6| ) is that the former appears near the exclusive boundary of the phase 



space (i.e. x — > 1,6 — > 0) while the latter occurs at a point x 1 = x c inside the physical 
region < x 1 < 1. The Sudakov singularities that arise at an exclusive boundary of the 
phase space are the most common and extensively studiedQ They are due to the loss of 
balance between the virtual contributions and the radiative tail of the real emission, which 
is strongly suppressed in this extreme kinematic regime. Sudakov singularities inside the 
physical region of the phase space have attracted less attention but are nonetheless quite 
common in jet physics. Some examples are the the C-parameter distribution [2I| in e + e 
annihilation and the jet shape [22] in hadron collisions. 

The general origin of Sudakov singularities at a critical point inside the physical region 
has recently been discussed 0. They arise whenever the observable in question has a non- 
smooth behaviour in some order of perturbation theory at that point. This can happen 
if the phase-space boundary for a certain number of partons lies inside that for a larger 
number, or if the observable itself is defined in a non-smooth way. Both mechanisms are 
responsible for the singular behaviour in the case of the isolated-photon cross section. The 
double- and single-logarithmic divergences in the square brackets of Eqs. ( f4.2| ), fl4.5|) are 
due to the first mechanism, while the remaining single-logarithmic terms in Eqs. ( |4.5|) 
are due to the second mechanism. We discuss these points in turn. 



"See, for instance, Refs. M, E0| and references therein. 
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At the LO the phase-space region r 7 < 1 is not accessible to the fragmentation compo- 
nent because there is a single parton, the triggered quark, inside the isolation cone. There- 
fore the coefficient function C^ L °)(x,r 7 ) in Eq. ( fOl ) has a step at the critical point r 7 = 1. 
As proved in Ref. J7|, this step- wise behaviour necessarily produces double- logarithmic sin- 
gularities at the next order in perturbation theory. When applied to our case, the general 
argument of Ref. is as follows. The NLO subprocess contains a term that is obtained by 
radiating a soft (and collinear) gluon from the LO subprocess. This term gives the following 
contribution to the NLO coefficient function 

C^ q NLO \x,r„6)= I' dzw{z)C[ LO \x^{z)) + ... , (4.9) 



where 1 — z is the energy fraction of the soft gluon (1 — z <C 1) and the dots stand for 
less singular terms when r 7 — > 1. In Eq. (|4.9|) , w(z) denotes the probability (in units of 
a 5 /27r) of emission of the soft gluon from the triggered quark; it is given by the usual 
'plus'-distribution in Eq. (O). The integration over z follows from energy conservation 



and relates the value of the parameter r 7 after the emission (i.e. on the left-hand side) to 
the corresponding value r 1 {z) before the emission. To explicitly evaluate Eq. ( f4.9|) , we have 
to specify how r 7 (z) depends on r 7 and z. Since r 7 is constrained by the fraction of the 
total energy inside the isolation cone (see Eq. ( |2.11| )), its value decreases (r 7 < r^ ut \z)) 
or increases (r 7 > r^ n \z)) according to whether the soft gluon is radiated outside or inside 
the cone. Neglecting less singular termsQ in the soft limit 1 - z< 1, we can thus write 

r (out)(^ ^ r7 + (i_2) ? (4.10) 
r^\ z ) ~ r 7 -(l-z) . (4.11) 



Using the explicit expression ( |3.3| ), ( |4.7| ) and inserting Eqs. ( |4.10[ ) and fl4.11| ) in the integral 
]9|), we respectively obtain: 

c (nlo, out) f Xjr s) = 2 8{l-x)C F f 1 dz In—*— ) 9(r 7 - z) + ... 

Jo VI — z 1 — zj + 

= + 6(1 -r 7 )5(l -x)C F ln 2 |r 7 - 1| +... , (4.12) 



Cf LOM \x ) r 11 5)=2 5{l-x)C F f dz (- In- ) 9(r 7 - 1 - (1 - z)) + . . . 

' JO \ 1 — Z 1 — Z J -(- 

= -0(r 7 - 1)5(1 -x)C F In 2 |r 7 - 1| +... , (4.13) 
in agreement with the double-logarithmic terms in Eqs. (|4.2|), 



Equations ( |4.9j ), ( f4.12| ), ( [4.13| ) explain the mechanism that produces the logarithmic 



divergences in the square brackets of Eqs. (|4.2j ), ( f4.5|) . The 'plus'-prescription in Eq. ( [4.7| 



arises from adding real and virtual soft-gluon radiation and, as a result of the cancellation 
of the soft singularities, leads to finite quantities whenever it acts on smooth functions 
of z. However, this is not the case of Eq. (|4.9| ), because C^ LO \x, r y (z)) is discontinuous 
at r 1 {z) = 1. The divergent behaviour at the critical point is thus due to the imperfect 
compensation || between real and virtual contributions, which occurs in the presence 



tThe actual size of the coefficient in front of the shift (f — z) on the right-hand side of Eqs. ( 4.10 ), ( 4.11 ) 
would affect only the single-logarithmic contributions. 
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of the LO step-like discontinuity at x 7 = x c . Because of the different kinematic recoil in 
Eqs. ( |4.10D and ( 4.11|) , at NLO the cross section has a double-sided singularity (Fig. [I]), that 



is, it diverges to +00 and —00 below and above the critical point x 1 = x c , respectively. 

The single-logarithmic term[] in Eq. Q4.5|) (outside the square bracket) and that in 
Eq. (|4.6|) have a different origin. They arise from the integration of the collinear spec- 
trum of the parton that is radiated by the triggered parton a at an angle #i 7 inside the 
isolation cone. They have indeed the form: 

P M ^ = P M Hmax + lnQ 2 /V) , (4.14) 

Jfi 2 /Q 2 tf l7 v J 

where P qa {x) is the relevant Altarelli-Parisi probability 

P qq (x) = C F ( , P qg {x) = C F l + {l ~ X? . (4.15) 

The lower limit of integration over comes from the factorization of the collinear sin- 
gularity at 6*i 7 = in the non-perturbative fragmentation function. The upper limit 6*^ ax 
comes from the kinematics of the process. Of course, #i 7 has to be smaller than the cone 
size 5. However, when r 7 — > 1 at fixed x and 5, the energy radiated inside the cone by 
the splitting process q — > a can violate the isolation constraint before the cone boundary 
is actually appproached by #i 7 . Therefore, 6^ ax ~ r 7 — 1 and Eq. ( |4.14 ) gives the single- 



logarithmic terms of Eqs. (|4.5| ) and ( |4.6| ). These terms are due to the non-smooth character 
of the isolation criterion, which enforces sharp boundaries on the energies and angles of the 
radiated partons. 



5 Outlook: higher orders, resummation 
and non-perturbative effects 

In this paper we have shown that the factorization of collinear singularities in cross sections 
for the production of isolated particles defined by the criterion ( [L.l|) is valid to any order 
in QCD perturbation theory. The non-perturbative component of the scattering process 
that is not power-suppressed at high transferred momentum Q is thus taken into account 
by the universal fragmentation functions D H / a (x,Q 2 ), whereas the isolation condition is 
consistently included in the short- distance subprocess. 

In the case of isolated photons produced in e + e~ annihilation, we have checked the 
factorization pattern by performing an explicit calculation at NLO in a^. This calculation 
shows that, although infrared and collinear safe, the short-distance component of the NLO 
cross section has still a divergent behaviour (Fig. [I]) when the photon energy fraction 
x 1 approaches a critical value x c = 1/(1 + £/,) that is located inside the physical region 
< x 7 < 1. As shown in Eqs. ( |4.2| ), (|4.3|), (|4.5|), (fL6|), the divergences are double- 



logarithmic. They are due to the loss of balance between real and virtual contributions 



*We recall that these single-logarithmic divergences appear also in the LO coefficient function C 7 of 
the direct component of the cross section. 
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that is enforced by the non-smooth character of the energy isolation criterion. In Sect. [| 
we have discussed in detail the physical mechanisms that produce these singularities of the 
Sudakov type. 



The same mechanisms leading to Eqs. (fO|), (fO|), (}4.5|), (|4.6|) will enhance the double- 
logarithmic divergences by further integer powers of In | r 7 — 1 1 in yet higher orders of 
perturbation theory. Since these divergences are unphysical 0, QCD calculations at any 
finite order in perturbation theory cannot give reliable phenomenological predictions for 
the isolated-photon cross section in the region around x-y = x c . 

This problem may be overcome by trying to avoid the phase-space region near the 
critical point. However, also in this case, some general theoretical understanding of the 
phenomenon is necessary to assess the extent of the dangerous region. The identification 
of the dangerous region in the perturbative calculation is even more difficult in the case of 
hadron collisions ||. Here the singularities appear at a critical point in the partonic cross 
section and are smeared by the convolution with the parton distributions of the colliding 
hadrons. 

A better way to deal with the problem is to identify and resum the soft-gluon divergences 
to all orders in perturbation theory. The resummation approach has been used successfully 
for the treatment of the Sudakov singularities near the exclusive phase-space boundary for 
many observables |19], pOfl . The same approach was advocated in Refs. [| for isolated 



photons and in Ref. J7[ for a general treatment of the Sudakov singularities inside the 
physical region. 

All-order resummation has already been carried out for a particular type of critical 
point inside the physical region. This type of critical point (Fig. H) regards cross sections 
that at LO have a stepwise behaviour at that point and at NLO show a divergence on 
a single side of the step. After resummation, the cross section is finite, continuous and 
differentiable at the critical point. Rather than suppressing the NLO divergence and thus 
mantaining the step-like behaviour, resummation leads to the suppression of the step. The 
general form of the resummed cross section is a smooth extrapolation from the region 
where the NLO divergence appears towards the critical point, joining smoothly with the 
finite value that the cross section has on the other side of the LO step. This characteristic 
structure was called a Sudakov shoulder [0]. 

The resummation of the double- sided singularities (Fig. [XJ) of the isolated-photon cross 
section is technically more complicated. We expect that it leads to a smooth 'jump' struc- 
ture obtained by smearing (convoluting) two Sudakov shoulders: a shoulder below and an 
inverted shoulder above the critical point x 7 = x c . As a result, the two sides of the LO 
step match at some intermediate value at the critical point in the all-order distribution, 
and the step is smoothed. The actual derivation of the resummed calculation requires more 
study, since the type of smearing (convolution) to be applied to the two shoulders strongly 
depends on the detailed isolation kinematics. Work on resummation is in progress, and the 
results will be reported elsewhere. 

Note that the perturbative divergences at x 7 = x c correspond to integrable singularities 
and therefore they could in principle be removed by non-perturbative smearing effects as 
expected, for instance, from hadronization. However, since the hadronization smearing 
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X 

c 

Figure 2: Divergent behaviour at a critical point with single-sided singularities: LO (solid), 
NLO (dotted). The dot-dashed line denotes the Sudakov shoulder obtained after all-order 
resummation. 

should cancel divergent terms proportional to some power of as(Q), this would require 
that the short- distance cross section contains non-perturbative contributions that are not 
power-suppressed at large Q. On the basis of our resummation argument, we do not 
anticipate the presence of these contributions. The resummation of soft-gluon effects to 
all orders of perturbation theory should be sufficient to render the isolated-photon cross 
section finite and smooth throughout the physical phase space. This suggests that the 
non-perturbative contributions that are not included in the fragmentation function are still 
power-suppressed. 
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